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Abstract 

Results of a numerical study concerning flow past a 70-deg blunted 
cone in hypersonic low-density flow environments are presented using the 
direct simulation Monte Carlo method. The flow conditions simulated are 
those that can be obtained in existing low-density hypersonic wind tunnels. 
Results indicate that a stable vortex forms in the near wake at and below a 
freestream Knudsen number (based on cone diameter) of 0.01 and the size 
of the vortex increases with decreasing Knudsen number. The base region of 
the flow remains in thermal nonequilibrium for all cases considered herein. 

Nomenclature 

A — base area of cone, ltd 2 /A 

C D = drag coefficient, s 2D/p 00 V^A 

C}{ = heat transfer coefficient, 2q/p 00 V^ ) 

d = base diameter 

D — drag 

Kn = Knudsen number, X/d 
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M = Mach number 

p — pressure 

R = radius 

Re — Reynolds number, pVd/p 

Re 2 = total Reynolds number, Re c<) (p 00 / p 0 ) 

s — distance along the body surface measured from stagnation 

point 

S = speed ratio 

T — thermodynamic temperature 

Tt = translational kinetic temperature 

u — axial velocity 

v = radial velocity 

V = velocity 

x = axial distance from stagnation point measured along symmetry 
axis 

{xjd) 0 — location of rear stagnation point 
( x/d)\ = location of wake stagnation point 
y = radial distance from symmetry axis 

A = size of the vortex, = (x/d) i — ( x/d) a 
A = mean free path 

(i — dynamic viscosity 

p — density 

Subscripts 

a — axial 

c — corner 

n — nose 

stag = forebody stagnation point values 

r = radial 

w — surface values 

o — stagnation chamber values 

oo = freest re am values 


Introduction 

Aeroassisted Space Transfer Vehicles (ASTV) will have three primary 
components: aerobrake, payload compartment, and main propulsion unit. 
These vehicles will maneuver within the rarefied and generally undefined 
outer atmosphere rather than just passing through it. Recent studies 1-4 
in the transitional flow regime of the aerothermal loads of various blunt 
bodies associated with current and projected ASTVs have focused mainly 
on the aerobrake. However, little is known about the wake structure of 
the aerobrake, which will interact with the payload compartment and main 
propulsion unit. The heating rate and aerodynamic forces that can result 
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from the interaction of the near wake with the payload compartment can 
influence significantly the size and shape of the afterbody configuration. 
There is a lack of experimental data for blunt body wake flows in the 
hypersonic rarefied flow regime. The structure of the near wake of an 
aerobrake in the continuum flow regime has been numerically studied 
recently, 5 using the Navier-Stokes equations. Although the forebody may 
be in the continuum flow regime, the expansion into the wake region can 
produce large local Ivnudsen numbers. Consequently, it becomes imperative 
that t he wake be investigated with an analysis capable of properly modeling 
the general features of low density- flows. 

Recent direct simulation Monte Carlo (DSMC) simulations 6 of rarefied 
hypersonic flows about spheres show that the calculated sphere drag agrees 
well with measured wind-tunnel values. These calculations also show that 
the w T ake is very rarefied, and that there are no vortices in the wake for 
freestream Knudsen numbers between 0.8 and 0.009. For similar flow 
conditions, Brewer’s 8 DSMC blunt body calculations also agree with these 
results. The Navier-Stokes calculations of Ref. 7 for similar flow conditions 
agree well with experimental drag but predict a steady vortex in the wake at 
a freestream Knudsen number value of 0.009. Such discrepancies highlight 
the importance of using a particle simulation method such as DSMC to 
establish a database for isolating the features of wake flows under conditions 
where rarefaction effects are significant. Furthermore, parallel studies using 
Navier-Stokes and DSMC simulations can establish the bounds for which 
continuum methods are appropriate. 

In the present paper, the DSMC method 9 ’ 10 has been used to simulate 
the low-density axisymmetric flow of nitrogen about a 70-deg half-angle, 
blunted cone with a base diameter of 5 cm (a wind-tunnel model of the 
Viking lander aeroshell 1 1 ). This model provides a realistic generic ASTV 
configuration for investigating the wake structure. 


Conditions for Calculations 

The flow conditions considered are attainable in the SR3 tunnel at 
Centre National de la Recherche Scientifique (CNRS 12 ), and are listed 
in Table 1. The model size selected is such that it can be tested in 
existing low-density hypersonic wind tunnels (Table 2). The surface 
temperature is assumed to be constant at 300 K, and the gas-surface 
interaction model is diffuse with full thermal accommodation. Because 
of the relatively low total temperature, the DSMC simulations model the 
nitrogen flow as a nonreacting gas with energy exchange between rotational 
and translational modes. The freestream parameters along with selected 
results are summarized in Table 2. 
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a Model diameter d = 4.95 cm, nose radius: R n = 1.25 cm, comer radius: R c = 0.0.036 cm 
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Results and Discussion 

Attention is focused on the flowficld structure of the forebody and the 
wake regions, but with major emphasis on the wake. In addition, surface 
quantities and drag resulting from low-density simulations are presented 
for the flow about a 70-deg blunted cone. 

Forebody Flow Structure 

Figure 1 shows the density contours for case 2. Near the stagnation 
point, a substantial density increase occurs for all three cases considered. 
This is characteristic of hypersonic flows about cold blunt surfaces. It has 
been observed from the density contours for the two most rarefied cases 
that the density rises gradually as the flow approaches the body, indicating 
the diffuse nature of the shock wave that is characteristic of highly rarefied 
flows. However, for case 3, where KtIoq = 0.001, the density contours 
show a rapid rise of density as the flow approaches the body, indicating the 
formation of a thick shock wave. In addition, the calculated results indicate 
(not shown) that the Mach 1 contour approaches the body surface at the 
aft corner, and the forebody is enveloped in a rather thick subsonic layer. 



Fig. 1 Nondimensional density contours (Case 2, A/oo = 19.7, A*Uoo = 0.0107, 

/ifoo = 2220). 
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Wake Structure 

The density in the near wake for all three cases is less than 40 % of 
the freestrearn value. The density even drops below 10 % of the freestream 
value in a region near the base, and this zone of very low density increases 
in size as the freestream Knudsen number increases. It has been also 
observed from these calculations that there is a substantial amount of 
thermal nonequilibrium in the entire wake region for cases 1 and 2. For 
case 3, most of the wake region is in thermal equilibrium. There is a small 
region near the base where the flow remains in thermal nonequilibrium. In 
addition to these thermal nonequilibrium effects, a substantial portion of 
the wake flow in the base region is subsonic, and the extent of this subsonic 
region increases with decreasing freestream Knudsen number. 

T he streamline plots show that the flow is attached to the surface, and 
there is no evidence of flow separation for case 1. This is expected because 
of the rarefied flow condition (Kn^ = 0.0317 and Re <*, = 768). However, 
for cases 2 and 3, the flow separates beyond the corner on the base of the 
body and a stable vortex develops in the base region of the wake (shown in 
Fig. 2 for case 2 only). It can also be inferred from these results that the 
separation point moves closer to the body corner with decreasing freestream 
Knudsen number. 

The wake stagnation point is generally defined for separated flows as 
the location of a point in the wake where the flow reattaches and the flow 
velocity is zero. For axisymmetric separated flow, this point lies on the 
symmetry axis of the wake where the axial velocity is zero. The effect of 
rarefaction on the location of the wake stagnation point is shown in Fig. 
3. The wake stagnation point moves toward the rear stagnation point with 
increasing freestream Knudsen number and finally coincides with the rear 
stagnation point for attached flow. Therefore, the size of the vortex (A) 
increases with decreasing rarefaction. The axial and radial components of 
translational temperature along the wake symmetry axis for cases 2 and 
3 are shown in Fig. 4. The translational thermal nonequilibrium exists 
all along the symmetry axis for case 2. For case 3, the flow along the 
wake symmetry axis is in translational equilibrium downstream of the wake 
stagnation point, and the translational temperature attains the highest 
value at the wake stagnation point (Fig. 4). 

Figure 5 presents the axial velocity variation in the radial direction at 
three different axial locations in the wake for case 3. These results show that 
the shear-layer thickness changes very gradually with distance downstream 
of the base, and that the leading-edge bow shock wave gradually becomes 
weak as the flow moves downstream in the wake. In addition, the lip 
expansion fan diffuses in the streamwise direction and there is no indication 
of any distinct wake compression wave. For the more rarefied cases, it is 
observed that there are no distinct wake flow features, as shown in Fig. 5 
for case 3. The lip expansion fan merges gradually with the shear layer and 
bow shock as the freestream Knudsen number increases. 
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Fig. 6 Effect, of rarefaction on the surface heat transfer. 


Surface Quantities 

The surface heat transfer coefficient is presented in Fig. 6 for all three 
cases. The results are shown as a function of the nondimensional distance 
{s/ R n ) along the surface measured from the forebody stagnation point. The 
expansion of the flow about the aft corner for case 3 causes the maximum 
heat transfer at that point, rather than at the stagnation point (Fig. 6). 
The present calculations also show that the surface pressure is not very 
sensitive to rarefaction since the forebody is enveloped with a fairly thick 
subsonic layer. The surface quantities on the base surface are very small 
compared with those on the forebody surface for all three cases. 

Conclusion 

Direct simulation Monte Carlo simulations of blunt body flows have 
been performed for flow conditions that can be achieved in existing low- 
density wind tunnels. For the range of the Knudsen numbers considered 
(0.03-0.001), the calculated results show that a stable vortex exists in the 
base region of the wake for freestream Knudsen numbers between 0.01 
and 0.001, and the size of the vortex increases with decreasing freestream 
Knudsen number. There is no formation of a lip separation shock or a 
distinct wake shock at these rarefied conditions. This is contrary to the 
behavior observed for blunt body flows at continuum hypersonic conditions. 
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